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Abstract 

The paper contains an application of the generalized lattice model to multicomponent systems 
with internal degrees of freedom. The short-range inter-atomic repulsions and smooth long-range 
parts of the inter-atomic potentials are considered separately by means of packing condition and 
in effective field approximation, respectively. The dependence of the inter-atomic potentials on 
the internal degrees of freedoms (such as atomic electric and/or magnetic momentum) taken into 
account. The Helmholtz free energy functional in the generalized lattice model is reduced to 
the Ginzburg-Landau-Cahn-Hilliard-like (GLCH) form. The connection between the inter-atomic 
potentials characteristics and the parameters of the GLCH-like functional is obtained. Equations 
for both equilibrium and quasi-equilibrium states in condensed systems are derived. It is shown 
that equilibrium distribution of the fast internal degrees of freedom by frozen space distribution of 
the components obeys to the Schrodinger-like equation. 
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I. INTRODUCTION 



The generalized lattice model (GLM) of muiticomponent condensed systems (such as 



account the following essential factors: 

1. The short-range inter-atomic repulsions. These repulsions are not identical for different 
pair of atoms, therefore it is impossible to take into account the repulsions by means 
of lattice introduction. 

2. The presence of the local fields due to the long-range parts of inter- atomic poten- 
tials. These fields have the essential influence on both equilibrium properties and 
non-equilibrium processes on the corresponding space scales. 

3. The Helmholtz free energy functional contains the well defined parameters that have 
connections with characteristics of the components and their interactions. 

The present paper contains further development of GLM. In addition to previous results 
it takes into account the following: 

1. Existence of the internal atomic degrees of freedom such as atomic electric and mag- 
netic moments. These degrees of freedom are responsible for the local magnetization 
and local electric polarization in the system. 

2. Existence of some non-equilibrium degrees of freedom due to colossal time of their 
relaxations. Real condensed systems duTQ . clS db rule, partially equilibrium systems. 

3. Existence of essential inhomogeneities due to both difference of the atomic sizes of the 
components and presence of some frozen degrees of freedoms. It is well known, the 
inhomogeneities scales can be arbitrary up to nanocluster sizes. 

II. SEPARATION OF SHORT-RANGE AND LONG-RANGE PARTS IN THE 
INTER-ATOMIC POTENTIALS 

Short-range and long-range parts of the inter-atomic potentials in condensed matter play 
essentially different roles. Short-range repulsions prevent the collapse of particles into multi- 
particle conglomerates and lead to some restriction on the local densities nj(r) of the com- 
ponents 



where u>i is the inverse value of the maximal local density of i-th component (i = 1 4- m, 
m is the number of the components in the system). The quantity Ui has dimensionality of 
volume and henceforth will called as the specific atomic volume of z-th component. As far as 
the quantity cjjnj(r) is the local volume fraction of z-th component at the point r, then we 
have the restrictions on the local densities of the components for all points r in the system: 





(1) 



m 




(2) 
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For condensed matter the inequality in this condition should be replaced on the equality: 

m 

rk(r) -1=0. (3) 

i=i 

This constrain (the packing condition) of the components local densities in condensed matter 
takes into account the short-range inter-atomic repulsions. It will be used in generalized 
lattice model (GLM) instead of the hypothesis on the lattice structure existence. The 
Helmholtz free energy minimization for equilibrium state of systems should be realized under 
the packing condition, which takes into account the short-range parts of the inter-atomic 
potentials in the system. Hence the inter-atomic potentials should be included into the free 
energy with cutting out their short-range parts: 

KM = { ^ (r) ' " I" 1 (4) 
n ' \ 0, otherwise, v ; 

where Wij{v) is "true" interaction potential between z-th and j-th components, are the 
cutting parameters, related to the specific atomic volumes of the components by the relations 

Oy ~ [(^) 1/3 + (^) 1/3 ] • (5) 

In addition of the packing condition Q, the numbers of the components atoms in the system 
should be fixed at the free energy minimization. These conditions have the following form: 

(V) 

The packing condition ([3]) and condition of numbers particles conservation ([H]) should be 
satisfied for any form of the configuration part of the free energy. 

III. HELMHOLTZ FREE ENERGY IN THE GENERALIZED LATTICE MODEL 

The Helmholtz free energy functional of a system with account of the atomic electric Dj (r) 
and magnetic moments Mj (r) of the components in presence of the external electric E (r) 
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and magnetic H (r) fields has the following form: 



(7) 



F = \ J2 J J K ^ r ~ ^>i{r)nj{r')drdr' 
E // - r >i(rK(r') (D«(r) ■ D,(r')) drdr* 
+5 2 // ^ " r ')^(r)^(r') (M,(r) • M,(r')) drdr' 
4 £ // ^ r - r ')^(r)^(r') (D,(r) ■ M,(r')) drdr' 

m „ m „ 

+ Z) / (E(r)-D. J (r))n i (r)rfr + ^ / (H(r) ■ M. 4 (r)) n,(r) rfr+ 

* =1 (V) i=1 (V) 

1-1 (V) 

where (r), Qjj (r), i?^- (r) and S^- (r) are the long-range parts of the relevant two-body 
inter-atomic potentials related to moments independent, electric dipole - electric dipole, 
magnetic dipole - magnetic dipole and magnetoelectric terms, respectively, the last term 
in this relation is the entropy term, T is the temperature in energetic units, n (r) is the 
summarized local density of the particles 

rn 

™(r)=$>,.(r) (8) 

i=i 

The equilibrium distributions of the components in the space are determined by mini- 
mum of the Helmholtz free energy functional with account of the conditions ([3]) and ([6]). 
But beyond these conditions, some connections between electric Dj(r) and magnetic M,(r) 
moments and the external fields should be included. Suppose that external fields influence 
on the orientations of related moments but do not influence on their magnitudes Di and Mj. 
Hence we have 

(Pi (r)) 2 - D\ = 0; (Mj (r)) 2 — M? = 0. (9) 

IV. THE LAGRANGE FUNCTIONAL 

For the minimization of the Helmholtz free energy (jTJ) at the conditions ([3]), ([6]), (Q, 
let us introduce the Lagrange functional £ depending on the local densities nj(r) of the 
components, their electric and magnetic moments Dj(r), Mj(r), the external fields E(r), 
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H(r), and the Lagrange multipliers ^(r), Aj(r), Vi(r), 



m 
i=l 



m(r) dr - Ni 



(V) 



m „ / \ „ / m 

E /^n?(r) [(Mj (r)) 2 — M 2 ] dr - j ^) 

i=l (y) (y) V i! 



-E J^nUr) [(D, (r)) 2 - A 2 ] * 

' 1 00 " (10) 



cjj rij(r) - 1 I dr, 



where F is the Helmholtz free energy defined by (JTj). 

The necessary condition of the Helmholtz free energy extremum is vanishing of the partial 
and functional derivatives of the Lagrange functional with respect to n,(r), Dj(r), Mj(r), 
^(r), Aj(r), ^i(r), Calculations of these derivatives lead to a system of nonlinear integral 
equations with kernels Kij (r — r') , Qij (r — r') , Rij (r — r') , Sij (r — r') . Unfortunately, at 
present there are no effective methods of such kind equations solutions with kernels of general 
form. But under the some conditions this system of equations can be reduced to a system of 
partial differential equations. Instead of the interatomic potentials (r), (r), (r), 

(r) this system of partial differential equations contains a set of integral characteristics 
of the potentials. 



V. REDUCTION OF THE GENERALIZED LATTICE MODEL TO GINZBURG- 
LANDAU-CAHN-HILLIARD-LIKE APPROXIMATION 

One of the most effective methods in statistical thermodynamics of condensed matter 



based on the phenomenological Ginzburg-Landau-Cahn-Hilliard (GLCH) models [ll|, [12 
Unfortunately, prognostic capabilities of the GLCH models are restricted by absence of the 
direct connections between model parameters and inter-atomic potentials. 

The GLM contains the inter-atomic potentials in the explicit form, but the mathematical 
structure of the free energy functional for GLM is more complicate than in GLCH model. 
By some additional assumptions, the GLM can be reduced to the GLCH-like model. 

There are at least three scales of the sizes in the system: 

1. atomic sizes a ; 

2. range of actions of long-range parts of the inter- atomic potentials ro; 

3. distances bo on which changes local compositions and/or local moments of the com- 
ponents in the system. 

Suppose these parameters obey the inequalities: 

a ^ r < b . (11) 

Let us expand the functions rij (r'), JDj (r'), Mj (r') in vicinity of point r in powers of 
r — r' up to second order terms. Using this expansion and the Green formula 

J u(r) Av(r) dr = - J (Vu (r) • Vv (r)) dr, (12) 

(V) (V) 
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that is valid if functions u (r), v (r) and their gradients Vw (r), Vi> (r) vanish on the boundary 
of the domain V, we obtain the following expression for Lagrange functional: 



1 m f f 

c = -uH / 4° ^ « ■ V ^ ( r )) 



E ( v \Pi (*) ^ ( r )] ■ v i D i w n i Ml) 



0=1 
3 



+itf £ (V [M? (r) (r)] • V [Mf (r) (r)]) 



a=l 



+24 2) £ (V [A° (r) n t (r)] • V [M° (r) nj (r)] ) W 



1 m 

4e 



a=l 
'(0) . ^(0) 



+ Qg' (D, (r) • D, (r)) + Rf (M, (r) ■ M, (r)) 



(V) 



+2S® (D.W-M.W) 



rii (r) (r) dr 



+E 



i=i 



(E (r).D* (r)) + (H (r)-M, (r)) + Tln 



(V) 



rii r 



n (r) 



E 

i=l 



/i, 



nj(r) dr - 



E 

i=l 



A,- (r) 



[(D, (r)) 2 - A 2 ] n 2 (r)> 



00 



- E / ^ [(M, (r)) 2 - M 2 ] ^ 2 (r) dr - y *(r) £ 

f=i ( V) (y) \i=0 



Ui nAr) - 1 \ dr, 



where 



K 



(p) 



Kij (r) |r| p dr, 



(p = 0, 2), 



00 



(13) 



(14) 



are the integral characteristics of the inter-atomic potentials; the parameters , Rf) , 



are defined similarly. 

Functional (Tl3|) is similar to Ginzburg-Landau and Cahn-Hilliard (GLCH) functionals, 
but in contrast to GLCH the integrand in (1131) is not a polynomial over order parameters 
and all the parameters have clear physical sense, due to their explicit connections (|14p with 
inter-atomic potentials. This functional is the sum of two terms 



C = £i + £ 



2' 



(15) 



m D 2 f M 2 f f 

A = E t** N i + -j- I X ^ n i ( r ) dr + ~y I ^ (r) n ^ (r) dr + I * {r)dr (16) 

i=1 [ (V) (V) J (V) 

is the part of the functional ({TBI , which does not depend in explicit form on functions Dj(r), 
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Mi(r), ni(r), and 



£2 = -12 E / K S } ( r ) • ( r )) 

ij=l (£) L 

3 

+Qif E ( V (r) m (r)] • V [Df (r) n, (r)]) 

3 
a=l 

+2Sif (V \D? (r) n, (r)] • V [Mf (r) n, (r)] ) U 



1 m 

^>E 

«J=1 



a=l 



(0) 



(V) 



+2S® (D.W-M.W) 



rii (r) rij (r) dr 



-E 



E 

i_1 (v) 



(E (r) • D, (r)) + (H (r) • M, (r)) + T In 



rii r 



n r 



n 2 (r) D 2 (r) + ^ n 2 (r) M 2 (r) + ¥(r) W , n«(r) \ dr 



(V) 



(17) 



is the second part, depending on n^r), Df(r), Mf (r) in explicit form. 

The term £ x does not give any contribution into derivatives of the functional (fT3|) with 
respect to rij(r), Df{r), M"(r). Therefore, the variational problem for the functional (|T3|) 
is equivalent to the variational problem for the reduced functional £2 with additional con- 
ditions ([3]), (El), dU). This reduced functional has a form 



(18) 



£ 2 (K(r)}) = / A(w s (r), Vu s (r)) dr, 

(v) 
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where u s (r) denotes all the functions n;(r), Df(r), Mf(r), and 



AKW. Vw a (r)) = -— J] 



< 2) (Vn 4 (r)-Vn,(r)) 



ij'=l 



E ( V [A Q (r)n,(r)] • V [^(rjn^r)]) 

a=l 
3 



+2Sf (V [A Q (r)n,(r)] • V [Mj*(rH(r)]) 



3 



0=1 



(19) 



4- + w ( D -M ■ D j( r » + < ( M <w ■ m jM) 



8=1 



+25^ (D,(r)-M,(r)) 
n 2 (r) 



rii[T)nj(T) 



{(E(r) ■ D,(r)) + (H(r) • M,(r))} n,(r) + ^ [A,(r)D?(r) + i*(r)M?(r 



+ E 

i=l 



T In 



ft r 



nir 



/i j - *(r) w< 



Thus, the equilibrium distributions of the components and their electric and magnetic 
moments obey the Lagrange-Euler system of equations for functional (ITS!) 



OA 



duJr) 



V 



dA 



d(Vu s (r)) 



0, 



(20) 



together with conditions (j3J), (jHJ), (EJ). 

Substitution (111) into (EEO) with u s (r) = m(r), u s (r) = Df (r), u a (r) = Mf (r) leads to 



S 



following equations: 

Kf) + Q<°> (D,(r) ■ D,(r)) + fl|°> (M,(r) ■ M,(r)) + sf) (D,(r) ■ M^r)) 



nj{r) 



m 3 



+ EE ^? ( VD °W • V [D;(rK(r)]) +4 2) (VM? (r) ■ V [M«(r)n,(r)]) 

[4 2) (VA°(r) ■ V [M«(r) Kj (r)]) + S« (VMf (r) • V [^(r)n,-(r)] )" 
(E(r) ■ Di(r)) + (H(r) • M,(r)) + n t (r) [A,(r)D 2 (r) + z/,:(r)M 2 (r 



+ 



+ 



+T In 



HAT 



nir 



(21) 



^ + t^* (r) 



1 m 



4° An,-(r) + £ qJ 2) (V • { Df (r) V (£>«(r) n,(r))}) 

+ E 4 2) ( v • ( M f « v ( M ;( r ) ^w) }) 

a 

+ E^ }< | ( V - {A ? W V(M?(r)n J -(r))}) + (V • {Mf (r) V (D?(r) n,(r)) }) 



a=l 



m 

£ {q§> D°( T ) + 4°) M/(r)} n,(r) r*,(r) + E a (r) n,(r) + A,(r)n 2 (r) Df (r 

1 m 

"6 E ( V ^« • V [W^WD + S if ( V ^W • V [M;(r) % (r)])] 

j'=i 

= " 6 V • E {^to V [^( r )%(r)] + 4 2) ^W V [Mf^n^)] } 



(22) 



£ |4 0) Mf(v) + i 4 } D?(r)| n,(r) n,-(r) + # a (r) n,(r) + i/ 4 (r)n? (r) M? (r 

i m r i 

j=i 

7 V ' E {^00 V [ M ;( r )^( r )] + Sg^(r) V [D;(r)n,-(r 



(23) 



6 ^ I 



respectively. 

The system of equations ( J2T1 l22l 123]) with conditions (EJ [6], [9]) describes space distri- 
butions of the components and local electric and magnetic moments for the case of full 
thermodynamic equilibrium in the system. 



A. Partial equilibrium in the system 

Complete thermodynamic equilibrium in condensed systems can be realized, as a rule, on 
astronomical scales of time only. Suppose that the local magnetization and the local electric 
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polarization in solid-state structures are the fast variables while the space distributions of the 
components are the slow variables. Then the magnetizations and the electric polarizations 
can be considered on appropriate time scales as equilibrium functions while the space distri- 
butions of the components are frozen variables. There are many examples of such systems. 
In particular, it could be considered a layered system such as multiferroic magnetoelectric 
composites [13l. Il4j]. For this situation the functions rii (r) are the fixed by the sample prepa- 
ration technology while the functions Dj (r) and Mj (r) obey the equations (J22J) and (l2"5j) 
with conditions Q). 

Let us consider as example a model of two-component system. Let the first component 
particles have the electric moment D and second component particles have magnetic moment 
M. Space distribution of these components is prescribed by their local densities rt\ (r) and 
n 2 (r), respectively. The functional (fTT|) for this case has the following form: 



£ 2 ({D (r)} , {M (r)}) = — £ / {q (2) (V [D a (r) n x (r)] ■ V [D a (r) n x (r)]) 

Q= V) 

+,R (2) (V [M a (r) ni (r)] • V [M a (r) m (r)]) 
+2S {2) (V [D a (r) m (r)] ■ V [M a (r) n x (r)])\dr 

Q(°) (D(r)) 2 (m(r)) 2 + (M(r)) 2 (n 2 (r)) 2 



(V) 



(24) 



+2S (0) (D(r) • M(r)) m(r)n 2 (r) 



+ 



dr 



(D(r) • E(r)) n a (r) + (M(r) ■ H(r)) n 2 (r) 



(V) 



A(r) 



n 2 (r) (D(r)) 2 - ^n 2 (r) (M(r)) 



dr, 



where Q« = Qg J , i2« = S« = S2\ (p = 0, 2). 

Equilibrium distributions of electric and magnetic moments at given densities ni(r) and 
n 2 (r) obey the equations 

f Q {0) nl(r)D a (r) + S (0) m(r)n 2 (r)M a (r) + E a {r)m{r) - A(r)fi?(r)D Q (r) 
~ {Q {2) (Vm(r) ■ V [D a (r)m(r)]) + ^ (Vm(r) • V [M Q (r)n 2 (r)])} 
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1. 



(25) 



= ~Q V • {<2 ( V(r)V [F(r) ni (r)] + ^V(r)V [M Q (r)n 2 (r)]}; 

R {0) n\{r)M a {r) + l S (0) n 1 (r)n 2 (r) J D a (r) + # a (r)n 2 (r) - z/(r)n 2 (r)M Q (r) 
-i {i? (2) (Vn 2 (r) ■ V [M a (r)n 2 (r)]) + S (2) (Vn 2 (r) ■ V [£> a (r)m(r)])} 

= —V ■ |i? (2) n 2 (r)V [M Q (i> 2 (r)l + S (2) n 2 (r)V [£> a (rWr)]}. 
6 

These equations can be slightly simplified. With account of the conditions ([9]) the complete 
system for equilibrium distribution of the electric and magnetic moments has the following 
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form: 



o 



--^ 2 >A - 
6 



u 



u 



A(r 

» + 

^K(r)) 2 - J D 2 n?(r) = 0; 

a=l 
3 

]>>«(r)) 2 -M 2 n 2 (r) = 0, 
where A is the Laplace operator, and 



r) + 
1 



_i 5 (2) A _ 5 (0) 

6 



-R®A - + „( r ) 
6 



v a (r) - E a (r) = 0; 
v a {r) -H a {r) = 0; 



t a (r) = D a (r)m(r) 
> Q (r) = M Q (r)n 2 (r). 

The first and second equations in (I26p can be written in the matrix form: 

'[-igWA - Q(°) + A(r)] [-fS^A - S^] 

[-1^(2) A - 5C°)] [-|i?( 2 )A - + u(t) 















l-l 





(26) 



(27) 



(28) 



In absence of the external fields this matrix equation is fully similar to the coupled Schrdinger 
equations for functions u a (r) and u a (r). The last two equations in ( 1261) are similar to the 
wave functions normalization; they make possible to find the unknown functions A(r) and 
z/(r) in just the same way as the normalization condition in quantum mechanics permits 
to find the spectrum. Thus, classical statistical thermodynamics of the quasi-equilibrium 
condensed systems with internal degrees of freedom can be based on the standard quantum 
mechanical methods. 



VI. CONCLUSION 

The paper contains a development of the generalized lattice model in following ways. 

1. The separation of the short-range inter-atomic repulsions and the smooth long-range 
parts of the inter-atomic potentials. Short-range repulsions take into account by means 
of packing condition, long-range parts take into account in effective field approxima- 
tion. 

2. Dependence of the inter-atomic potentials on the atomic internal degrees of freedoms 
(such as atomic electric and/or magnetic momentum) taken into account. 

3. The Helmholtz free energy functional in the generalized lattice model is reduced to the 
Ginzburg-Landau-Cahn-Hilliard-like form. The connection between the inter-atomic 
potentials characteristics and the parameters of the GLCH-like functional is obtained. 

4. The equations for both equilibrium and quasi-equilibrium states in condensed systems 
are derived. It is shown that equilibrium distribution of the fast internal degrees of 
freedom by frozen space distribution of the components obeys to the Schrodinger-like 
equation. 
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One of the most interesting applications of the GLM is research of layered structures with 



alternation of magnetic and ferroelectric layers [13|, [14J. In particular, this approach permits 



to take into account lattices misfits on the interphases boundaries without any auxiliary 
assumptions; 

to impart the physical interpretation to the phenomenological models like GLCH mod- 
els; 

to find the ways for prognosis of the layered structures. 
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